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(Dated: June 24, 2005) 

We present a search for electroweak production of single top quarks in the s-channel and t- 
channel using neural networks for signal-background separation. We have analyzed 230 pb -1 of 
data collected with the D0 detector at the Fermilab Tevatron Collider at a center-of-mass energy 
of 1.96 TeV and find no evidence for a single top quark signal. The resulting 95% confidence level 
upper limits on the single top quark production cross sections are 6.4 pb in the s-channel and 5.0 pb 
in the t-channel. 



PACS numbers: 14.65.Ha; 12.15.Ji; 13.85. Qk 
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INTRODUCTION 

Top quark physics provides fundamental knowledge 
of the strong and electroweak sectors of the standard 
model and offers discovery potential for physics beyond 
the standard model. The top quark was discovered in 
1995 at the Fermilab Tevatron Collider in tt events pro- 
duced through the strong interaction 0|. The standard 
model predicts that proton-antiproton collisions should 
also produce single top quarks through the electroweak 
interaction. Studying single top quark production will 
provide direct measurements of the CKM matrix element 
|Vtb| and top quark polarization, and will probe possible 
new physics in the top quark sector 0, 0] . 

There are two main modes of single top quark pro- 
duction as shown in Fig. ^ the s-channel (tb) process 
pp — > tb+X and the t-channel (tqb) process pp — > tqb+X. 
The production cross sections have been calculated at 
next-to-leadingorder (NLO) in the strong coupling con- 
stant 0, H UM @, yieldmg0.88 ± 0.14 pb for the s- 
channel and 1.98±0.30 pb [4j,|5( for the t-channel, assum- 
ing a top quark mass of mt = 175 GeV. Both the D0 and 
CDF collaborations have previously performed searches 
for single top quark production |9j, yj| . Recently, CDF 
performed a search using 160 pb" 1 of data and obtained 
upper limits of 13.6 pb (s-channel), 10.1 pb (t-channel), 
and 17.8 pb (s+t combined) at 95% C.L. 




FIG. 1: Representative Feynman diagrams for electroweak 
top quark production at the Tevatron Collider. This figure 
shows (a) the s-channel and (b) the t-channel. 

In this Letter, we present a new search for electroweak 
production of single top quarks with the D0 detector. 
Our search focuses on the final state in which the top 
quark decays into a b quark and a W boson, and the 
W boson then decays into an electron or a muon, and 
a neutrino. This gives rise to an event signature with a 
high transverse momentum (j>t) lepton with significant 
missing transverse energy (fir) from the neutrino. In ad- 
dition, the s-channel has two 6-quark jets, whereas the 
t-channel typically has one 6-quark jet together with a 
light-quark jet. The 6-quark jet in the t-channel, as rep- 
resented in Fig. mb), is usually emitted in the forward 
direction with low pr and is often undetected Hi. The 
main backgrounds in this analysis are the W boson pro- 
duction in association with jets (W+jets), top quark pair 
(tt), and multijet production. We use neural networks to 



separate the signals from the backgrounds. In the ab- 
sence of any significant evidence for signal, we build a 
binned likelihood from the neural network outputs to set 
upper limits on the single top quark production cross 
sections. 



THE D0 DETECTOR 

The D0 detector consists of a central tracking system, 
liquid- argon/uranium calorimeters, and an iron toroid 
muon spectrometer [l2T ]. The central tracking system 
covers the detector pseudorapidity region |?7det| < 
2.5. It includes a silicon microstrip tracker and a scintil- 
lating fiber tracker, both located within a 2 T solenoidal 
magnet. The calorimeters consist of a central barrel in 
the region |?7det| < 1-1, and two end caps extending the 
coverage to |??dct| ~ 4.2. The muon system outside the 
calorimeter consists of a layer of tracking detectors and 
scintillation counters before 1.8 T toroids, followed by 
two similar layers after the toroids. Luminosity is mea- 
sured using plastic scintillator arrays located in front of 
the end calorimeters. 



DATA SET AND EVENT SELECTION 

The data were recorded between August 2002 and 
March 2004 using a trigger that required an electromag- 
netic energy cluster and a jet in the calorimeter for the 
electron channel, or a muon and a jet for the muon chan- 
nel. The integrated luminosity is 226 pb -1 for the elec- 
tron channel and 229 pb -1 for the muon channel. 

In the electron channel, we require exactly one iso- 
lated electron Q with pt > 15 GeV and |?/dct| < LI. 
In the muon channel, events are selected by requiring 
exactly one isolated muon Q with p T > 15 GeV and 
|?7dct| < 2.0. For both channels, events are also re- 
quired to have Jfir > 15 GeV. Events must have from 
two to four jets with the leading jet pr > 25 GeV and 
|?7det| < 2.5, and all other jets having pr > 15 GeV and 
|*7det| < 3.4. Jets are defined using a cone algorithm with 
radius 1Z = 0.5 15] . Misreconstructed events are rejected 
by requiring that the direction of fir is not aligned or 
anti-aligned in azimuth with the lepton or the jets. 

The fraction of signal-like events is further enhanced 
through the selection of 6-quark jets that are identified by 
reconstructing displaced vertices from long-lived parti- 
cles. A displaced vertex is selected by requiring the trans- 
verse decay-length significance, L xy /aL xy , to be greater 
than seven, where L xy is the decay length and ctl is the 
uncertainty on L xyi calculated from the error matrices of 
the tracks and the primary vertex. A jet is considered 6- 
tagged by this algorithm if a displaced vertex lies within 
a cone of radius 1Z = 0.5 around the jet axis 
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For both s-channel and t-channel searches, we separate 
the data into independent analysis sets based on final- 
state lepton flavor (electron or muon) and 6-tag multi- 
plicity. To take advantage of the different final state 
topologies, we separate single-tagged (=1 tag) events 
from double-tagged (>2 tags) events. In the t-channel 
search, we additionally require that one of the jets is not 
b tagged. 

ACCEPTANCES AND YIELDS 

We estimate the kinematic and geometric acceptances 
for s-channel and t-channel single top quark production 
using the COMPHEP matrix element event generator 0] 
with m t = 175 GeV. The factorization scales are ml for 
the s-channel samples and (mt/2) 2 for the t-channel sam- 
ples. For the s-channel (t-channel) search, the t-channel 
(s-channel) is considered as background. 

We use both Monte Carlo events and data to esti- 
mate the background yields. The VF+jets and diboson 
(WW and WZ) backgrounds are estimated using events 
generated with ALPGEN [J^. The diboson background 
yields are normalized to NLO cross sections computed 
with MCFM 0]. The total VF+jets yield is normalized 
to the yield in data corrected for the presence of multi- 
jet, tb and dibosons before requiring a 6-tagged jet. The 
fraction of heavy-flavor (Wbb) events is obtained using 
the ratio of the NLO cross sections for W+jets and Wbb, 
as described in Ref . [2(| • This normalization to data also 
accounts for smaller contributions such as Z+jets events 
in which one of the lcptons from the Z boson decay is 
not reconstructed. 

The ti background, consisting of the leptonic decay 
modes of the W boson from the top quark decay (£±jets 
and dilepton), is estimated using samples generated with 
ALPGEN, normalized to the cross section: a(tt) = 6.7 ± 
1.2 pb |21|. where the uncertainty on the top quark mass 
is incorporated into the cross section uncertainty. 

The parton-level samples are then processed with 
pythia 22] for hadronization, particle decays, and mod- 
eling of the underlying event. The generated events are 
processed through a GEANT-based |23( simulation of the 
D0 detector. The resulting lepton and jet energies are 
smeared to reproduce the resolutions observed in data. 

The background from jets misidcntificd as electrons or 
jets resulting in isolated muons is estimated using mul- 
tijet data samples that pass all event selection cuts, but 
fail the requirement on muon isolation or electron qual- 
ity [LH . This background is normalized using a data sam- 
ple dominated by multijet events, selected by requiring 
$ T < 15 GeV. 

The overall acceptances, including trigger and selection 
efficiencies, for signal events with at least one 6-tagged jet 
are (2. 7 ±0.2)% in the s-channel and (1.9 ±0.2)% in the 
t-channel. The acceptance is calculated as the fraction of 



events that pass the selection over all possible single top 
quark decays, including all leptonic and hadronic decays 
of the W boson. Estimates for signal and background 
yields and the observed numbers of events after selection 
are shown in Table [I] 

TABLE I: Estimates for signal and background yields and 
the numbers of observed events in data after event selection 
for the electron and muon, single-tagged and double-tagged 
analysis sets combined. The VF+jets yields include the di- 
boson backgrounds. The total background for the s-channel 
(t-channel) search includes the tqb (tb) yield. The quoted yield 
uncertainties include systematic uncertainties taking into ac- 
count correlations between the different analysis channels and 
samples. 



Source 


s-channel search 


t-channel search 


tb 


5.5 ± 1.2 


4.7 ±1.0 


tqb 


8.6 ±1.9 


8.5 ±1.9 


W+jets 


169.1 ± 19.2 


163.9 ±17.8 


ti 


78.3 ±17.6 


75.9 ± 17.0 


Multijet 


31.4 ±3.3 


31.3 ±3.2 


Total background 
Observed events 


287.4 ±31.4 
283 


275.8 ±31.5 
271 



NEURAL NETWORKS ANALYSIS 

After event selection, several variables are combined in 
neural networks to discriminate the single top quark sig- 
nals from the backgrounds. The networks are composed 
of three layers of nodes: input, hidden, and output. For 
training and testing, we use the mlpfit [24| package. 
Testing and training event sets are created from simu- 
lated signal and background samples. We use a tech- 
nique called early stopping to determine the maxi- 
mum number of epochs for training which prevents over- 
training. Each network is then tuned by choosing the 
optimal number of hidden nodes. From studies based 
on optimizing the expected upper limits on the single 
top quark production cross sections, we find that the 
s-channel and t-channel searches each require only two 
networks, corresponding to the dominant backgrounds: 
Wbb and ti ->■ ^±jets. 

The list of discriminating variables has been chosen 
based on an analysis of Feynman diagrams of signals and 
backgrounds [13] and on a study of single top quark pro- 
duction at NLO 0, 0] ■ The input variables to each net- 
work are selected from this list by training with differ- 
ent combinations of variables and choosing the combina- 
tion that produces the minimum testing error and largest 
signal-background separation. Table [D] shows the vari- 
ables used for each signal-background pair. These vari- 
ables fall into three categories: individual-object kine- 
matics, global-event kinematics, and angular correla- 
tions. 
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TABLE II: Input variables for each neural network signal-background pair. Variable descriptions can be found in the text. 







Signal-Background Pairs 








tb 




tqb 




Variable 


Description 




Wbb 


tt 


Wbb 


tt 




Individual object kinematics 












^ / 1 ,-,4- 1 \ 

PTU etl taggedJ 


Transverse momentum of the leading tagged jet 




V 


V 


V 


— 


p T (jetl untaggcd ) 


Transverse momentum of the leading untagged jet 




— 


— 


V 


V 


PT(jet2 untagged ) 


Transverse momentum of the second untagged jet 




— 


— 


— 


V 


p T (jetl nonbcst ) 


transverse momentum or the leading nonbest jet 




V 




— 


— 


Pr(jet2 nonbcst ) 


Transverse momentum of the second nonbest jet 

Global event kinematics 




V 


V 


— 


— 




Transverse mass of the two leading jets 




V 


— 


— 


— 


PT(J ctl . J et ^J 


Transverse momentum of the two leading jets 






— 


V 


— 


A'T ( 1 1 1 Q+D^ 

ivi ^aiijcts ) 


Invariant mass of all jets 




V 


V 


V 


V 


Jrljp ^aiijcts J 


Sum of the transverse energies of all jets 




— 


— 


V 


— 


M ^alljets - jetl taggcd j 


Invariant mass of all jets excluding the leading tagged jet 










V 


_h l^aiijcts — JCTi taggcd ; 


Sum of the energies of all jets excluding the leading tagged jet 




— 


V 


— 


V 


ff T (alljets - jetl tagged ) 


Sum of the transverse energies of all jets excluding the leading tagged jet 




— 


— 


— 


V 


p T (alljets - jetl tagged ) 


Transverse momentum of all jets excluding the leading tagged jet 




— 




— 


V 


AJ (alljets - jct bcst ) 


Invariant mass of all jets excluding the best jet 




— 


v/ 


— 


— 


.H (alljets - jet bcst j 


Sum of the energies of all jets excluding the best jet 




— 


V 


— 


— 


-Ht (.alljets - jet bcst ) 


Sum of the transverse energies of all jets excluding the best jet 




— 


V 


— 


— 




Invariant mass of the reconstructed top quark using the leading tagged jet 




V 


V 


V 


V 


M(W,jet be8t ) 


Invariant mass of the reconstructed top quark using the best jet 




V 


— 


— 


— 


VS 


Invariant mass of all final state objects 

Angular variables 






— 






A7?.(jetl,jet2) 


Angular separation between the leading two jets 




V 




V 




^U^J-untaggodJ X V£ 


Pscudorapidity of the leading untagged jet X lepton charge 








V 


V 


cos(£, Q^Xz)top bost 


Top quark spin correlation in the optimal basis for the s-channel |2(1 , reconstructing 


V 










the top quark with the best jet 












cos(£, jetl untagged )top taggod 


Top quark spin correlation in the optimal basis for the t-channel L'tl , reconstructing 






V 




the top quark with the leading tagged jet 












cos(alljets, jetl tagged ) a iij cts 


Cosine of the angle between the leading tagged jet and the alljets system 


in the 






V 


V 


alljets rest frame 












cos(alljets,jct nonbcst ) alljcts 


Cosine of the angle between the leading non-best jet and the alljets system 


in the 




V 







alljets rest frame 



Since the input variables do not depend on the lepton 
flavor, the electron and muon analyses utilize the same 
variables. However, owing to different lepton resolutions 
and pseudorapidity ranges, we construct separate net- 
works for them. Therefore, four neural networks are used 
for the signal-background pairs (tb-Wbb, tb-tt, tqb-Wbb, 
tqb-tt) for each of the electron and muon channels. 

Figure|21shows distributions of four representative vari- 
ables. We reconstruct the final state top quark from the 
reconstructed W boson and a jet as follows. The W bo- 
son is reconstructed from the isolated lepton and the 
missing transverse energy. The z-component of the neu- 
trino momentum {p v z ) is calculated using a W boson mass 
constraint, choosing the solution with smaller \p v z \ from 
the two possible solutions |2^. In the s-channel analysis, 
the top quark is reconstructed from the W boson and the 
"best" jet 0. The best jet is defined as the jet in each 
event for which the invariant mass of the reconstructed 
W boson and the jet system is closest to m t = 175 GeV. 
In the i-channel analysis, the top quark is reconstructed 
from the W boson and the leading 6-tagged jet. Using 
these two methods we are able to correctly identify the 



6-quark jet from the top quark decay in about 90% of the 
signal events. 

Figure |21 shows the outputs of the neural networks for 
the data and the expected backgrounds, as well as sig- 
nals for the electron and muon channels combined. The 
neural network output in mlpfit is around one for sig- 
nal events and around zero for background events, but it 
is not constrained to the interval [0, 1]. The tt networks 
separate signal and tt backgrounds efficiently. The Wbb 
networks are less efficient for the IU+jets backgrounds 
because the event kinematics are similar between signal 
and background. 

SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are evaluated for the Monte 
Carlo signal and background samples, separately for the 
electron and muon channels and for each 6-tag multiplic- 
ity. The most important sources of systematic uncer- 
tainty are listed in Table ITTT1 The systematic uncertainty 
on the shapes of the distributions is also taken into ac- 
count for the contributions from 6-tag modeling, jet en- 
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FIG. 2: Comparison of signal, background, and data for the electron and muon channels combined, requiring at least one tag, 
for four representative neural network input variables. Shown are (a) the transverse momentum of the leading untagged jet, (b) 
the invariant mass of the reconstructed top quark using the leading tagged jet, (c) the invariant mass of the final state system, 
and (d) the pseudorapidity of the leading untagged jet multiplied by the charge of the lepton. Signals are multiplied by ten. 



ergy calibration, jet identification, and trigger modeling. 
In order to evaluate the total uncertainty, we consider 
all sources of systematic uncertainties for all samples and 
their correlations. The total uncertainty on the signal ac- 
ceptance for single-tagged events is 13% for the s-channel 
and 15% for the i-channel, and for double-tagged events 
it is 24% for the s-channel and 28% for the i-channel. 
The total uncertainty on the background yield is 10% for 
the single-tagged samples and 26% for the double-tagged 
samples. 



CROSS SECTION LIMITS 

The observed data are consistent with the background 
predictions for all eight analysis channels. We therefore 
set upper limits on the single top quark production cross 
section separately in the s-channel and i-channel searches 
using a Bayesian approach |2^. In each search, two- 
dimensional histograms are constructed from the Wbb 
vs. tt neural network outputs. A likelihood is built from 
these histograms for signal, background, and data, as a 
product over all channels (electron and muon, single and 



TABLE III: Range of systematic uncertainty values for the 
various Monte Carlo signal and background samples in the 
different analysis channels. 



Source of 


Uncertainty 


systematic uncertainty 


range (%) 


Signal and background acceptance 




6-tag modeling 


5-20 


jet energy calibration 


1 - 15 


trigger modeling 


2-7 


jet fragmentation 


5-7 


jet identification 


1 - 13 


lepton identification 


4 


Background normalization 




theory cross sections 


2-18 


IT+jets flavor composition 


5-16 


Luminosity 


6.5 



double tags) and bins. We assume a Poisson distribu- 
tion for the observed number of events in each bin and 
a flat prior probability for the signal cross section. The 
prior for the combined signal acceptance and background 
yields is a multivariate Gaussian with uncertainties and 
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tb-Wbb NN output 




tqb-tf NN output 




tqb-Wbb NN output 



FIG. 3: Comparison of signal, background, and data for the neural network outputs, for the electron and muon channels 
combined, requiring at least one tag. This figure shows (a) the tb-ti filter, (b) the tqb-ti filter, (c) the tb-Wbb filter, and (d) the 
tqb-Wbb filter. Signals are multiplied by ten. 



correlations described by a covariance matrix. Finally, 
we compute the posterior probability density as a func- 
tion of the production cross section. 

The Bayesian posterior probability densities are shown 
in Fig. 01 for both the s-channel and ^-channel searches. 
The corresponding upper limits at the 95% C.L. are 
6.4 pb in the s-channel and 5.0 pb in the i-channcl. 
The sensitivity of these measurements is given by the 
expected upper limits obtained by setting the observed 
number of events to the background prediction in each 
bin. The expected upper limits are 4.5 pb in the s- 
channel search and 5.8 pb in the ^-channel search. 



represent significant improvements over previously pub- 
lished results 0, El EI due to the larger data set as 
well as the use of a multivariate analysis technique to- 
gether with shape information from the resulting output 
distributions. They approach the region of sensitivity 
for models of physics beyond the standard model, such 
as a fourth quark-generation scenario or flavor-changing 
neutral-currents 3] . 



CONCLUSIONS 



No evidence is found for electroweak production of sin- 
gle top quarks in 230 pb -1 of data collected with the D0 
detector at y/s = 1.96 TeV. The data consist of events 
in the electron and muon final states with at least one 
6-tagged jet. We build binned likelihoods from the out- 
put of neural networks to set upper limits at the 95% 
C.L. The measured s-channel limit is 6.4 pb and the 
measured i-channel limit is 5.0 pb. These upper limits 
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FIG. 4: The Bayesian posterior probability density as a func- 
tion of the single top quark cross section for the s-channel and 
t-channel searches. 
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